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Abstract 

We study non-equilibrium noise in the transmission current through barri- 
ers in 1-D Luttinger hquids and in the tunneUng current between edges of 
fractional quantum Hall liquids. The distribution of tunneling events through 
narrow barriers can be described by a Coulomb gas lying in the time axis along 
a Keldysh (or non-equilibrium) contour. We show that the charges tend to re- 
organize as a dipole gas, which we use to describe the tunneling statistics. The 
dipole-gas picture allows us to have a unified description of the low frequency 
shot noise and the high frequency "Josephson" noise. The correlation between 
the charges within a dipole (intra-dipole) contributes to the high-frequency 
"Josephson" noise, which has an algebraic singularity at u) = e*V/h, whereas 
the correlations between dipoles (inter-dipole) are responsible for the low- 
frequency noise. We show that an independent or non-interacting dipole ap- 
proximation gives a Poisson distribution for the locations of the dipole centers 
of mass, which gives a flat noise spectrum at low-frequencies and corresponds 
to uncorrelated shot noise. Including inter-dipole interactions gives an ad- 
ditional 1/t^ correlation between the tunneling events that results in a |u;| 



singularity in the noise spectrum. We present a diagrammatic technique to 
calculate the correlations in perturbation theory, and show that contributions 
from terms of order higher than the dipole-dipole interaction should only af- 
fect the strength of the |a;| singularity, but its form should remain ~ |a;| to 
all orders in perturbation theory. A counting argument also suggests that the 
leading algebraic singularity at ujj should be oc |a; — wjp^"^ to all orders in 
perturbation theory. 

PACS: 72.10.Bg, 73.20.Dx, 73.40.Gk, 73.50.Fq, 73.50.Td 
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I. INTRODUCTION 



The noise spectrum in a two-terminal conductor in the absence of an apphed voltage is 
proportional to the conductance and to the temperature. This result was found experimen- 
tally by Johnson in 1927 ]l|| and explained theoretically by Nyquist in 1928 0. Such relation 
between equilibrium noise and conductance can be seen as a consequence of the fluctuation- 
dissipation theorem. The noise in the presence of transport (non-equilibrium noise) can also 
be related to transport coefficients for non- interacting systems , but now these transport 
coefficients, in the most general case, cannot be determined from conductance measurements 
alone. For interacting systems one should expect an even richer behavior, as new features 
in the noise should appear as a consequence of correlations due to interactions. In general 
the shape of the noise spectrum is determined by the dynamical properties of the system, 
which contain information about the excited states. Thus the noise spectrum is a powerful 
probe which allows us to study dynamics of strongly correlated systems. 

Interacting electron systems at 1-D form strongly correlated states -Luttinger liquids- 
whose properties are well understood. However, it has been very difficult to realize 1-D 
Luttinger liquids in experiments. This is because even a small amount of impurities cause 
the localization of the 1-D electrons and destroy the Luttinger liquids. Recently it was 
realized that another strongly correlated 1-D state - Chiral Luttinger liquid - exists on the 
edges of fractional quantum Hall (FQH) liquids. Due to its chirality {i.e. all excitations 
move in the same direction) and the lack of back scattering, a chiral Luttinger liquid cannot 
be localized by impurities. Thus it is possible to realize, in practice, extended 1-D systems 
through FQH states. Recently an IBM group studied the tunneling between two edges of 
filling fraction 1/3 FQH states. They found that the tunneling conductance has a power law 
dependence on temperature which is a characteristic property of (chiral) Luttinger liquids 
0-^. Their finding is consistent with the theoretical prediction a oc for the u = 1/3 FQH 
state 0. In this paper we will study the noise spectrum in the tunneling current between 
(chiral) Luttinger liquids. The noise spectrum carries rich information about dynamical 
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properties of (chiral) Luttinger liquids, which will help us identify such strongly correlated 
states in experiments. 

Recent studies of noise in non-interacting systems reveal that the noise spectrum con- 
tain features that come from the statistics of the tunneling particles [|T0|. These statistics- 
dependent features are not contained in the DC conductance. For Luttinger liquids, the 
tunneling particles sometimes carry fractional statistics and fractional charges. It is then 
very interesting to study the noise spectrum for tunneling between (chiral) Luttinger liq- 
uids, especially those features that come from the strongly correlated properties of (chiral) 
Luttinger liquids (such as fractional statistics and fractional charges). 

Two kinds of noise may appear in tunneling at a finite voltage V, the shot noise and 
the "Josephson" noise. The shot noise can be understood from a classical picture in which 
the average tunneling current is viewed as a result of many tunneling events. A tunneling 
event represents a single particle (which can be an electron or a charged quasiparticle) 
that tunnels through the junction. The spectrum of the shot noise is determined by the 
correlations between tunneling events. In this paper we always assume that the tunneling 
time is much shorter than the average spacing between two tunneling events. Under this 
approximation, we will ignore the retardation and model the tunneling by an instantaneous 
tunneling operator Tipl^ipn + H.c, which transfers particles between two reservoirs. The 
"Josephson" noise is related to the fact that the two systems connected by the junction have 
different chemical potentials. The quantum interference between wave functions on the two 
sides of the junction may cause a singularity at frequency u = e*V/hm noise spectrum (here 
e* is the charge of the tunneling particle). Such features near the "Josephson" frequency 
ujj = e*V/h are called "Josephson" noise. In this paper we will develop a language for 
non-equilibrium noise in 1-D Luttinger liquids which covers both the shot noise and the 
"Josephson" noise. 

We start with the Keldysh formalism, in which the tunneling events are described by 
a Coulomb gas of charges on a Keldysh contour. Under certain conditions the charges at 
different branches of the contour pair into dipoles (in this case the Coulomb gas is said to 
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be in the dipole phase). The dipoles correspond to the tunnehng events in the shot-noise 
picture. The non-interacting dipole approximation leads to a Poisson distribution for the 
separation of dipoles, which results in a white noise {i.e. a frequency independent noise) at 
low frequencies. However, for a finite voltage across the junction, we find that the dipoles 
have a non-zero dipole moment which leads to a long range interaction between dipoles. 
The dipole-dipole interaction gives rise to a non-trivial distribution of the tunneling events 
which induces a \uj\ singularity in the low frequency noise spectrum. The dipoles have 
finite size and the intra-dipole structures are found to be responsible for the high frequency 
"Josephson" noise, which appear as an algebraic singularity of the form \uj — cjjp^"^ in the 
noise spectrum within perturbation theory. 

The full expression for the singularity at zero frequency in the noise spectrum due to the 
dipole-dipole interaction is found to be 



where It is the average tunneling current and g contains information on the interactions in 
the Luttinger liquid (or filling fraction of the FQH states, in the case of chiral Luttinger 
liquids) . Because of the nonlinear dependence of It on uj , the strength of the singularity 
in the noise spectrum at zero frequency will also have a non-linear dependence on ujj ( 



g = 1, where one recovers the \uj\ singularity that appears to order in the transmission 
coefficient D [|ll|. The correlations in the case of non- interacting electrons come from the 
Pauli principle, which enters very simply in the formulation used in this paper through the 
language of bosonization. 

The paper is organized as the following. In section II we will review the bosonization 
scheme for 1-D fermionic systems. In section III we calculate the noise perturbatively. In 
section IV we use the non-equilibrium (Keldysh) scattering operator as a means to obtain 
a joint probability distribution for tunneling events. The tunneling events can be mapped 
into charges of a Coulomb gas, which tend to reorganize as a dipole gas. A non-interacting 




(1) 




The particular case of non-interacting electrons can be obtained with 
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dipole approximation leads to uncorrelated noise. Dipole-dipole interactions and correla- 
tions will be discussed in section V, which lead to an \uj\ singularity in the low frequency 
noise spectrum. In section VI a diagrammatic technique is presented that accounts for the 
correlations in a systematic way. We show the existence of the \uj\ singularity at zero fre- 
quency to all orders in perturbation theory. A counting argument also suggests that the 
leading singularity at Uj should remain of the form \uj — ujI"^^'^ to all orders in perturbation 
theory. 



II. TUNNELING IN 1-D LUTTINGER LIQUIDS 

In this paper we will study the effect of particle interactions in the noise spectrum of 
a 1-D conductor. The results for 1-D systems of interacting particles, or 1-D Luttinger 
Liquids, can be directly used to study noise in the tunneling current between two edge 
channels in the Fractional Quantum Hall (FQH) regime. Figure 1 displays the geometries 
we are considering here. Fig. la shows a 1-D channel connected to two reservoirs, with 
a weak link or tunneling barrier in the middle of the channel. Figs, lb and Ic show two 
configurations in which we can observe tunneling between edge channels. The configurations 
can be accessed experimentally using metallic gates placed on top of the 2-D electron gas. 
Applying a negative gate voltage depletes the electron concentration underneath the gate, 
causing the two branches of edge states to get closer, and thus enhancing the tunneling 
between the channels. Because in this configuration both edges form the boundary of the 
same QH liquid, there can be either electron or quasiparticle (carrying fractional charge) 
tunneling. By applying a sufficiently large gate voltage, one can obtain the situation in Fig. 
Ic, where the edges form the boundaries of two disconnected QH liquids, and thus only 
electrons can tunnel from one edge to the other. 

Both the interacting 1-D systems and the FQH edge states are best described in the 



bosonized language fT^JT^ . In the case of the interacting 1-D system, the electron operator 
can be written as: 



6 



n odd 



where the and 9 fields satisfy the equal-time commutation relations 

TT 

x), e{t, y)] = -i-sgn(x - y) . (3) 

The canonical momenta associated with (j) and 6 are then = ^dxO and He = ■^dxcj) 
respectively. The dynamics of and 9 are described by the Hamiltonian density 

'H = ^[9{dxct>f + -^{dxef] , (4) 



where the effect of interactions enters through g ||T2[. For repulsive interactions g < 1, 
whereas for attractive interactions g > 1. For non-interacting electrons (7 = 1. The electron 
propagator has a power law decay envelope, with the long-range behavior dominated by 
(V'"^(t,0)^/'(0,0)) oc t~^9+9-^)n. 

The presence of a weak link or a potential barrier in the channel gives an additional term 
in the Hamiltonian which can be expressed in terms of the bosonic fields and 9 []^. For a 
potential barrier located at a; = the perturbation is oc il)'^{t,x = 0)ilj{t,x = 0), which can 
be written (keeping only the most relevant term) as 

n,nt = r6ix)e'''^'^'^ + H.c. . (5) 

For a weak link one can also show that the perturbation is 

n^nt = T6{x)e'''f'^''''^+H.c. . (6) 

Using a rescaling = 2y^0 and 9 = -^6, the Lagrangian densities for the small barrier 
and weak link problems are respectively 

£ = -^[{dt9)^ - {djf] - T6{x) e'^^em ^ ^ ^ /^^ 

OTT 

and 

^ = ^[{dt^f - {d.Af] - e^^^^*'°^ + H.c. (8) 

OTT 



Now, for the FQH edge states, we can write the right and left moving electron and 
quasiparticle operators as '^ji^L{x,t) =: e^^v^"^^'^*^^'*) : , where g is related to the FQH bulk 
state. For example, for a Laughlin state with filling fraction u = 1/m we have g = m for 
electrons and g = 1/m for quasiparticles carrying fractional charge e/m. The (pR^L fields 
satisfy the equal-time commutation relations 

, 0R,L(t,2/)] = ±«7r sgn(x - y) . (9) 

The dynamics of is described by 

J^R,L = ^ 5x0R,L i±dt - vd^)(pR^L , (10) 

where v is the velocity of edge excitations (which we will set to 1). The same algebraic decay 
of the electron operator occurs in the edge states of the FQH effect, where we have a chiral 
Luttinger liquid with the exponent g directly related to the bulk state (for a review see Ref. 

Hi)- 

The tunneling between left and right moving branches can be written as Hfun = r\E'|^\E'ij + 
H.c. We can write, in terms of = 0/? + 0^, the following total Lagrangian density: 

C = ^m^Y - (9,.0)2] - T6{x) e^v^'^(*'°) + H.c. , (11) 

with satisfying [0(t, x), dt4>(t, y)] = iniS^x — y). 

The Lagrangian for in Eq. ([Tl| ) is exactly the same as the one for 9 in Eq. (|^) and, 
with g ^ 1/g, the same as the one for in Eq. (|^). It is this Lagrangian in Eq. ([TT|) that 
will be the basis of our work. A voltage difference between the two reservoirs at the ends 
of the 1-D channel, or between the edges of the QH liquid, can be easily introduced in the 
model by letting F — >■ re"""^"*, where ujq = ujj = e*V/h, with e* = e for electron tunneling 
and e* = e/m for quasiparticle tunneling. 

Notice that in order to obtain the coupling term we assume that we only have contribu- 
tions from X = for the tunneling operators. This is the case when the width of the barrier 
is narrow. Also, if the barrier is narrow, the time spent in the tunneling is small compared 
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to the spacing between tunneling events. Indeed, in this case we can speak of tunnehng 
events that occur at rather well defined time coordinates. 

Using this language, the average tunneling current through a barrier in a one dimensional 
channel and between edge states in the FQH regime was calculated The current 

has a nonlinear dependence on the applied voltage between the terminals, with the power 
dependence on the voltage intimately related to the exponent g in the electron propagator. 
For the case of tunneling through a single barrier in a 1-D channel, or non-resonant tunneling 
between FQH edge states, one finds that It ~ V"^^"^ at zero temperature. In this paper we 
will study the noise in this current, starting with a perturbative calculation and then moving 
to a formalism that grasps non-perturbative contributions. 



III. PERTURBATIVE CALCULATION 

We can show that the tunneling current operator is Itit) = jit) = ■je*re*v^'^*^*'°^ + H.c. . 
For example, in the case of tunneling between edges (such as in Figs, lb & Ic) we simply 
use that It = —^[Nl,II] = ■^[N^,H] (where Nr^i are the total charge operators on the 
R, L edges) and the commutation relations to obtain the expression for It. Similarly, we can 
find the same for the case of a 1-D interacting system. The noise spectrum can be obtained 
by calculating two-point correlations involving the operator It{t). 

Notice that, as the problem under consideration is intrinsically non-equilibrium, one 
should use the Keldysh (or non-equilibrium) formalism in computing expectation values 



of operators. This is the case here, where if we treat the coupling term perturbatively and 
introduce an adiabatic turning on and off of the interaction, the state at t = — oo differs 
from the one at t = oo; the charge transfer in one direction due to the applied voltage 
clearly makes the two states at ±00 different, as the total charge in one edge branch (or 
reservoir) decreases whereas in the other the total charge increases. This problem could, in 
principle, be circumvented by including another term in the Hamiltonian that would close 
the circuit and bring the charges that tunneled through the barrier back to the reservoirs 
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(a "battery"). Such a way of thought is relevant to clarify the distinction between the 
equilibrium and non-equilibrium formalism, and how they can be connected in principle. 
However, in practice, adding the restoring charge coupling in the Hamiltonian only would 
make the problem more cumbersome and poorly defined, which makes the non-equilibrium 
formalism a natural choice. 

For perturbative calculations of zeroth and first order, however, there is no difference 
between the results for expectation values obtained with either the equilibrium or the non- 
equilibrium formalism. This is the case in the calculation of the current-current correlation, 
where the lowest order contribution is the zeroth order: 

The non-zero contributions come from the terms that, when applied to |0), transfer zero 
total charge, so we can write 

(j(t)j(O)) = e*^|rp (e'''^''\0\e'^^^'^e-'^'f'^^^\0) +e''^'''{0\e''^'''e-'^^^^ 

The (f) field correlation is (O|0(i)0(O)|O) — —2\n{5 + it), where 5 is an ultraviolet cut-off 
scale. The current-current correlation is then given by 

«tM(o)> = e-|rp^^. (14) 

which displays clearly oscillations at frequency f — ujq/2t: — e*V/h. This implies that the 
noise spectrum will also display structure at this frequency. The noise spectrum is calculated 
from the current- current correlation: 

/oo 
rf^e-*({j(t),j(0)}) 
'OO 

= e*^|r|^[c+(a;o + a;) c_(a;o + a;) c+(a;o - u;) c_(a;o - u;)] , (15) 

where 

/•oo p-i^P Ott , ,^ 
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The c-i-(c<j) will appear again in the next section, where we shall obtain their finite temper- 
ature version. The noise spectrum to order |rp is then given by 

= e*It[\l- u/ujol'^-' + |1 + u/uol^'-' ] , (17) 
where we used the perturbative result to order |rp for the tunneling current It = 

From the expression for S{uj) above we can deduce some features of the noise to order 
|rp. First notice that for u <^ ujq we obtain S{uj) ~ 2e*It, the classical shot noise result, 
independent of g. Notice also the singularities at uj = ±uJo. In the particular case of non- 
interacting electrons {g = 1) we have S{lj) = 2e*It for < uq, and S{uj) = 2e*It\Lj\/uJo 
for |ti;| > uq, which agrees, to lowest order in the transmission coefficient D (lowest order 
in |rp), with previous results for the noise spectrum of non-interacting electrons [|ll|]. To 
get the term in we need to go beyond this zeroth order perturbation theory, as we will 
do later in the paper. The sharp edge of the noise spectrum at uj = e*V/h for g = 1 finds 
its origin in the Pauli principle, which is the sole factor responsible for correlations in the 
non- interacting case [|1^]. In our model, particle statistics enter automatically in the way we 
construct the electron/quasiparticle operator from the boson fields and their commutation 
relations. 

In the following sections we shall see how the low-frequency noise spectrum is modified 
once we go beyond this perturbative calculation. 

IV. THE JOINT PROBABILITY DISTRIBUTION 

As we have previously mentioned, when the tunneling barrier is narrow so that the time 
the charge spends in the tunneling process is small compared to the times between two 
consecutive tunnelings, one can speak of well defined tunneling events at certain specific 
times. In this section we will find a joint probability distribution for the times for these 
tunneling events. 
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The term e*''"^ in the Hamiltonian (where we use 7 = y^) transfers charge from one edge 
branch to the other (say, in the case of the geometry of Fig. lb & Ic, from the R to the L 
edge branch). The term e~'^'^'^ does the converse (from L to R). We will map the problem 
to a Coulomb gas in a 1-D space, associating a charge + to the term e^'^'^ and a charge — 
to e~^'^'^. Let Z = {0\Sc{—oo, —oo)\0), where Sc{— 00, —00) is the scattering operator in the 
contour from t = — (X)tot = oo, and back to t = — 00 (the Keldysh formalism contour). In 
terms of the usual scattering operator S, we can write: 

Z = (0|5(-oo, 00) ,5(00, -oo)|0) 

= (0|5t(oo, -00) S{oo, -oo)|0) . (18) 

In this form the contributions from the forward (t = —00 — >■ 00) and return (t = 00 — > —00) 
branches are easily identified in terms of the more commonly used (equilibrium) scattering 
operators. Clearly, since S is unitary, Z — 1. Now let's expand Z perturbatively in F. We 
will use the scripts t and b to denote the top (or forward) and bottom (or return) branches, 
and + and — to denote whether the inserted operator is e*''"^ (+) or e~'^'^'^ (— ). Q+t- wiU 
denote the number of times that e*''"^ or e""*''"^ appear in the top and bottom contours (see 
Fig. 2). With this notation, we can expand the scattering operator as: 

Q*+,Q*- W+- J i=i j=i i=i j=i 

and 

5'(— 00, 00) = S^{oo, —00) 

Q+,Q'L '^=1 '=1 '==1 '=1 

where T stands for reverse time ordering. Notice that in the operator S'c(— 00, —00) = 
5'(— 00,00) 5'(oo,— 00) the T ordering occurs to the left of the T ordering, so that we 
replace both by a Tc ordering operator such that times in the top branch are ordered in- 
creasingly, times in the bottom branch are ordered decreasingly, and times in the bottom 
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branch are always ordered after the ones in the top branch (see Fig. 2). Using Tc we can 
write SJ—oo, — oo) as 



E 

r+,Qt,Q''+,Q'L 



' QifgLJ^W f n n k n < n 



i=l j=l k=l 1=1 



i=l 

In order to calculate the bracket 



E 

i=i 



Q5 



fc=i 



(o|r,(e*^[ ^(t*+)-E.=i ^(tr)+E.=i ^(tfe+)-Ei="i i)|o) 



(20) 



we use 



(0|Tc(e*'''^^*^ eV0(«'))|o) = (.-ii'ioiTcmMt'm) 



(O|T,(0(tO0(t2))|O) 



(21) 



and the contour-ordered two-point correlation 

— 21n((5 + i\ti — t2\) , both ti and t2 in the top branch 

—2 ln(5 — i|ti — t2|) , both ti and ^2 in the bottom branch 

—2 ln(5 — ^(ti — t2)) , ti in the top and t2 in the bottom branch 

—2 ln((5 + i{ti — ^2)) , ^1 in the bottom and ^2 in the top branch 



The bracket in Eq. (^) contains the contributions from all pairs of charges, which 
interact via a two body potential that is determined by the ordered two-point correlation. 
The phase terms due to Uq (e"""^"* for a + charge, and e*"^"* for a — charge) correspond to an 
underlying background, which tends to polarize the gas, leaving (in the case of positive uo, 
for example) more -|- charges than — ones in the top branch, and more — charges than + 
ones in the bottom branch. An illustrative picture of the unbalance created by the applied 
voltage V (non-equilibrium) is shown in Fig. 3. One can think of V as an "electric field" 
that polarizes the Coulomb gas, leaving an unbalance of + and — charges in the t and b 
contours, which gives rise to a net current in one direction or the other (excess of +(— ) 
charges, or R L {L R) tunneling), depending on the sign of V. 
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The expression for Z obtained as an expansion in F is exact so far. Also, the map 
into a Coulomb gas model is now complete. An expansion similar to the one we present 
here appears in the study of dissipative quantum mechanics models in a periodic potential 



I5| , p!6[] . There the charges are grouped in terms of the so called "soujourns" and "blips". 
We find the idea of keeping the + and — charges more intuitive, as is the idea of having the 
non-equilibrium voltage be thought of as a "field" that polarizes the gas and changes the 
densities within the t and h contours. This language, as we will show, makes it easier for us 
to go beyond the independent blip approximation, and study correlations. 

We will now focus in showing how the expression for Z can be used to define a joint 
probability of tunneling events. In the limit of a narrow barrier, as we pointed out previously, 
one can speak of rather well defined tunneling times or tunneling events. In this limit we 
can interpret the times that enter in the perturbative expansion of Z as the times for real 
tunneling events, and the sums and integrations as the means of including all tunneling 
histories in a partition function. Notice that it is very important that we understand that 
this interpretation has a meaning only when the tunneling barrier is narrow. 

Also notice that only the tunneling times in the forward or top branch can have a physical 
interpretation as a tunneling of a real charge (we only observe increasing times, with the 
return branch being just a mathematical tool). The correct joint probability distribution of 
tunneling events should be obtained by integrating out all t^^'s. This is a difficult task, and 
we shall appeal to a more intuitive picture that will allow us to sort out the most important 
contributions. This more intuitive picture can be extracted from the Coulomb gas model 
depicted in Figure 4. 

The first step we take is to recast the sum in terms of dipole configurations, as opposed 
to a sum of charge configurations. The dipole is determined by a center of mass coordinate 
terra and a dipole strength p. There are four types of dipoles, as shown in Fig. 5. The type 
of dipole depends on which branches the + and — charges are located at. We call a t dipole 
one in which both charges are in the top branch. A h dipole is one where the charges are 
in the bottom branch. In a c+ the + charge is on the top and the — on the bottom. In a 
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c_ the converse is true, the — is on the top and the + on the bottom. This distinction is 
important, as we will see it later. 

For a given charge configuration labeled by {Q\_,Q^_,Q^,Q^_} we associate a dipole 
configuration {rit, rib, n^, where the n's are, respectively, the number of t, b, c+, and c_ 
dipoles. The n's and Q's are related by: 

QX = nt + n+ ^nt + n- 

Q+ ^ nb + n_ = rife + n+ . 



Rewriting Z in terms of the n's instead of the Q's becomes a simple combinatoric task, 
which gives: 



E 



.j)Qi+Qt (i)Q++Q- , . . ^ 

t ) ^ ^ :^:nt A^nt inb inb 



nt,rn„n+,n- 



X 















\ "- J 


\ "- J 



ntWn+\nJ. x INTEGRAL 



y ^ — ^ — ^^ 1 1 X INTEGRAL (22) 

where the INTEGRAL term contains the interactions between the charges integrated over all 
positions. The first approximation we will make is what we will call the "independent dipole" 
approximation. The attraction between opposite charges tends to bind them together, and, 
if the fugacity of the gas (measured by |rp) is small, we can to lowest order neglect the 
interaction between dipoles. The only interactions entering in the calculation of Z are the 
intra-dipole interactions. The INTEGRAL term in the dipole approximation can be factored 
as a product of the contributions of individual dipoles. 

INTEGRAL = 6"" C++ c!r (23) 

where 
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e 



— oo 



One can check that t + 6 = c+ + c_ , so that summing over all Ut and in Eq. (^) can be 
shown to yield: 

Z = ^ (E£^(E£_r:. (25) 

Let us now interpret this expression. As we mentioned above, only events occurring in 
the forward or top branches can be observed. Therefore, the occurrence of a dipole of the 
c+ type implies a tunneling event in one direction occurring at the vicinity of the center 
of mass coordinate of the dipole. Conversely, a dipole c_ implies a tunneling event in 
the opposite direction. The statistical distribution of these center of mass coordinates of 
dipoles appears in the noise. The uncertainty of the location of the charges comprising the 
dipole with respect to the dipole center of mass also contributes to the noise; this intra- 
dipole contribution, however, is already partly taken care of in the first order perturbative 
calculation of noise, which can be seen to be nothing but the correlation between the position 
of the two charge components of a dipole. The intra-dipole noise is in the high-frequency 
range, centered at u = ujq = e*V/h. The contribution to the noise that we obtain with the 
Z in Eq.(^) is in the low-frequency range {u <^ uq), where the positions of the charges and 
dipole centers are not distinguished. The reason why we summed over the dipoles of type 
t and b is that they do not contribute to the noise beyond the intra-dipole order. These 
types of dipole correspond to tunneling in one direction shortly followed by tunneling in the 
opposite direction, which contribute to noise in the time scale of the dipole size, included in 
the intra-dipole contribution. 

With the interpretation above in hand, we can use Eq. (pS]) to argue that, in the dipole 
approximation, the tunneling events in either direction are independent, with a distribution 
that is Poisson-like with two parameters: |rpc+ and |rpc_. The probability of tunneling in 
one direction in an infinitesimal time At is P+ = |rpc+ At, the probability of tunneling in 
the opposite direction is P_ = |rpc_ At, and the probability of no tunneling event in this 
time is 1 - (P+ + P_). 
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This two-parameter Poisson distribution can be used to reproduce the results obtained 
for the tunnehng current to first order in perturbation theory. The tunnehng current is 
simply given by It = e*|rp(c+ — c_), i.e., the net rate of tunneling in one direction. To 
obtain an expression for It in terms of V we need to evaluate c+ and c_ : 

c+K) = c4-^o) = y_ dpj^—-^ = 2^^^ . (26) 



We can obtain the finite temperature results for c± by a conformal transformation ||T7 



e 



■OO 



dPTT -T^ e*"3 / dp . ^ , , (27) 



which gives 



c4^o) = C-(-c^o) = 2(7rT)2^-^5(^ + ^,^?-|^) cosh(^) (l + tanh(^)) , (28) 

where B is the Beta function. Using these expressions for c± we obtain 

/, = Ae*\mnTr<^-'Big + ^,g-^) sinh(^) , (29) 

which is the same expression found by first order perturbation theory in Ref . . For T = 0, 
in particular, we find that It ~ e*|rpV"^^^^. 

We now turn to the noise properties derived from this dipole approximation. Because 
the distribution in this approximation is Poisson like (and therefore uncorrelated) , we should 
expect the noise to have a flat frequency dependence, i.e., white noise. We are left with 
the problem of determining the amplitude of the noise. For this purpose, we will follow 
a calculation similar to one presented by Landauer |T^. Let (j^)A/ be the component 
of the noise power spectrum that falls in the frequency interval A/. Let also 5'(/) = 
Jodtj{t) e-^"^*'*, where 6 is a time interval. These quantities are related by: 

0")a/ = lim A/ . (30) 

The charge transferred in a small interval of time r is ±e* (with probabilities |r pc-tr), or 
0. We can write j{t) = J2njo{t ~ '^'t) In, with g„ = ±1, 0. Here jo is a narrow current pulse 
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that fits a slot of time r (tlie widtli of the pulse should determine a cut-off frequency above 
which the spectrum is no longer fiat). We can then write 

Sif) = / dte-'-'J2Mt -nT) = J2 e~^""^gn / du e-^^^Uu) . (31) 

Jo „ J —nr 

The last integral can be approximated by the total charge that tunnels (e*), since the pulse 
is narrow compared to r. We then have S{f) = e* J2n^~^'^^^Qn, and thus 

\Sif)\' = e*' e-'"""e^""'"g„g„' • (32) 

n,n' 

The uncorrected tunneling implies that {qnqn') = {q)"^ + {{q"^) — {<i)'^)5n,n' ■ After summing 
over n and n' we obtain that |S'(/)p = e*^A^((g^) — (q')^), where = B/r is the number of 
time slots. Now (g) = |rp(c+ — c_) r and (g^) = |rp(c+ + cJ) r, and for small tunneling 
times compared to the time between tunneling (g) ^ 1, so that we can neglect (g)^ and 
obtain 

(/)A/ = 2e*'|r|2(c+ + c_)A/. (33) 

We can connect the white noise amplitude to the tunneling current using Eqs. (|28|) and 
(p9|), and obtain 

(j2)a/ = 2e*h coth(^)A/ . (34) 

If we write It = GV = Guo/e* and take the ujq ^ limit, we obtain = 4TGA/, 

which is nothing but the Johnson-Nyquist equilibrium (V^ = 0) result. The non-equilibrium 
white noise can then be cast in a simple relation to the equilibrium Johnson-Nyquist noise, 
which is 

{f)Af = C^)cothC^){jrif. (35) 

The expression above for T — > gives (j^)a/ = 2e*/j A/, which is the classical expression for 
shot noise. Quantum corrections to the shot noise only come to order |r|^ and higher, and 
thus do not appear in the independent dipole approximation (order |rp). Also notice that 



the expression connecting equilibrium and non-equilibrium noise ^^cotli(^^) is universal 
independent of g and thus independent of interactions. This is consistent with the fact 
that the independent dipole approximation is a lowest order perturbative result, so that the 
assumptions necessary for the fluctuation-dissipation theorem are satisfied. 

The dipole approximation therefore captures the uncorrelated part of the noise. In the 
next section we shall see how correlations come about. 



V. BEYOND THE INDEPENDENT DIPOLE APPROXIMATION 

In this section we shall improve the dipole approximation. We have seen that the location 
of the centers of mass of two dipoles is uncorrelated in the approximation of the preceding 
section. In order to observe correlations one must include in the model the interactions 
between distinct dipoles. This is the next order correction to the INTEGRAL term in Eq. 



Consider two dipoles as show in Fig. 6. We take them, for the sake of illustration, to be 
both of the c+ type. The INTEGRAL term can be written for this case as: 



/ dtidt2dpidp2 



-tUJOPl 



-IUJ0P2 
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X 



[5 + i{t2 + P2/2 -ti- pi/2)fs [5 - t{t2 - P2/2 - ti + pi/2)] 



2g 



-ILJ0P2 



dtidt2dpidp2- 



which after the p\ and p2 integration yields 



1 + 2(7 



(t2-ti)^ 



dt\dt2 — 2g- 



c'c' 



This can be generalized to any two types of dipole to 



C4.C+ / dtidt2 e 



(36) 



[5 + i{t2 - P2/2 - ti - Pi/2)]23 [5 - i{t2 + P2/2 - ti + pi/2)Y9 ■ 

For dipole separations that are large compared to dipole sizes ([^2 — ^i| ^ 1^21), we can 
expand the expression in the integrand to obtain 

P1P2 



(37) 



(38) 
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did2 / dtidt2 e ^ (*2-'i)' , (39) 

where di^2 can be any of t, 6, c+ or c_, and rf'^^ 2 stands for the derivative of di^2 with respect to 
ujq. Using a similar argument to the one we used to obtain the finite temperature expression 
for c±{ujq)^ we can obtain the finite temperature version of the dipole-dipole interaction by 
simply substituting t2 — ti by sinh[7rT(t2 — ti)]/ ij^T) and using the T 7^ results for c±{uq). 
Nevertheless, we will just concentrate for the rest of the paper on the T = problem. 

From Eq. (^) we read that the dipoles interact through a potential. This dipole- 



dipole interaction gives rise to a non-trivial distribution of tunneling events, which show up 
in the noise spectrum as a cusp at zero frequency. Before proceeding to obtain the explicit 
form, including the strength of the singularity, we must understand when this picture that 
the charges can be assembled in pairs starts to break down, and correlations not contained 
in this dipole picture become important. 

The assumption we made in order to obtain correlations as in Eq. ( ^9]) was simply that 
the dipole sizes were small compared to the separation between dipoles. The mean dipole 
separation is related to the average current /j, and is given by e*//j. The dipole size can be 
taken to be the d' /d in Eq. (|39[) , since it is this term that enters in the interaction between 
the dipoles, and thus measures the distance between the + and — charges that form the 
dipole (Notice that, because the charges in the Coulomb gas are ±1, the distance between 
the -|- and — charges equals the dipole strength). The expressions for t and h depend on the 
cut-off scale 5, whereas c± are finite as we take 5-^0 (we can show that c+ + c_ = t + 6, 
and that the divergences in t and 6, which are purely imaginary, cancel each other). We 
then have that t' /t and h' /h must both scale as 5, and c'^/c± = {2g — 1)uJq^ (using Eq. ( p8[ ) 
and setting T — > 0). Therefore the dipole approximation is good as long as cug ^ <^ / It, or 
It « {e*yh)V. 

In the case of tunneling between edge states, this is the limit of a small tunneling current 
as compared to the Hall current. In the case of the 1-D channel, this limit corresponds to 
a small tunneling current as compared to the current for the non-interacting case. Because 
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of the non-linear I — V characteristic of the tunnehng current in 1-D Luttinger hquids 
{If oc V^^~^), the cases g > I and g < I are quite distinct. For g > I the dipole phase 
exists for small applied voltages V, whereas for g < 1 the dipole phase exists for large V. 
Now, one can still use the results of the dipole phase to study the noise in the case of (? > 1 
and large V, and the case of g < 1 and small V, by resorting to the duality g \/g that 
connects the two configurations shown in Figs, lb & Ic. The idea is that as one increases 
the applied voltage between the R and L edges in the configuration shown in Fig. Ic, the 
tunneling current It increases asymptotically, tending to the Hall current. Deviations from 
the Hall current correspond to "defects", or tunneling in the direction perpendicular to the 
Hall current, which is exactly the direction of tunneling shown in Fig. lb. Similarly, one 
can go from the situation in Fig. lb to the one in Fig. Ic by decreasing the applied voltage 
between the R and L edges. The bottom line is that, by wisely choosing which current 
direction to focus on, one can most often place the problem in the dipole limit for either one 
configuration or its dual with g ^ 1/g. The regime in which the dipole picture fails is then 
at the crossover between the two configurations, where the gas will be in a plasma phase. 

Now that we understand when the approximation is valid, let us look at its consequence 
in the noise spectrum. At zero temperature we only have either one of c+ or c„ dipole types, 
depending on the sign of ojq (see Eq. (p6D ). For concreteness, let us take ujq > 0, so that 
c+ dipoles survive. Since t' /t^h' /h ~ 5, the main correlations come from the interactions 
between c+ dipoles (for voltages small compared to 1/5), so that for large times the density- 
density correlation for c+ dipoles (which equals the current- current correlation) is given 
by 

(p+(t) P^mc-{p^? ~^'%;^^\ (40) 

which gives a noise spectrum S{uj) = 4:ng(2g — 1)^ (^-^^ e'^^^^^/Xc, where Ac is a short time 
scale cut-off (of order ujq^) for the 1/t^ correlation. The leading singularity at low- frequency 
is then 

S,,ng,{u;) = Angi2g-lfi^f \uj\ . (41) 

LlJq 
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Since It oc V'^^~^, the strength of the singularity has a non-hnear dependence V"^(3~^) on the 
apphed voltage. 

For the particular case of non- interacting electrons {g = 1) one can write ^ = where 
D is the transmission coefficient. The noise spectrum singularity is then SgingX^) — ^~^\^\: 



which recovers the result of Ref. [|TT|. The effects of correlation due to the Pauli principle 
enter automatically in our formulation of the problem through the bosonization. 

To finish this section, let us consider the case of equilibrium noise within the interacting 
dipole approximation. For g > 1 the tunneling current vanishes for V = 0. In the case of 
g < I, however, we have to invoke the dual picture {g — > 1/g) in order to use the dipole 
language. In any case the \u! \ singularity due to the dipole-dipole interaction vanishes for V = 
0. The reason can be viewed very simply: the non-equilibrium voltage was responsible for the 
polarization of the dipole gas, and the dipole-dipole interaction gave the correlation. At 
equilibrium, the average dipole strength vanishes, and the interactions in this case must come 
from induced dipole, or "Van der Waal's" attraction, which for our log potentials goes as 
We can show that the low-frequency behavior of the noise no longer has the |ci;| singularity, 
but has leading contributions from u"^ and \uj\^. The leading singularity is then oc \uj\^. The 
contributions calculated above are only the inter-dipole correlations. We should also account 
for the intra-dipole correlations, because for V = the singularity at the "Josephson" 
frequency falls to u = 0. We already calculated the intra-dipole contribution to the noise 
spectrum perturbatively in section III, which for zero applied voltage is S{uj) oc |c(jp^^^. At 
low frequencies, the intra-dipole noise will be dominant for g < 2, while the inter-dipole 
noise will be dominant for g > 2. Notice that, in the equilibrium case, an expansion for 
Z like the one in Eq. (|I9|) could be carried out with only one branch, since in equilibrium 
there is no need for the two branches of the Keldysh contour. 
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VI. DIAGRAMMATIC TECHNIQUE 



The dipole gas picture we used to expand Z can be justified in more formal manner. 
In this section we shall present a systematic way to expand Z diagrammatically, which 



is used in one-dimensional dissipative quantum mechanics models |19]. In this expansion 
we can identify the terms we included in the dipole picture. The expansion is the formal 
support for the more intuitive and physical picture of the dipole gas. We will first present an 
introduction to the diagrammatic expansion, followed by the calculation for the equilibrium 
case and implications for the non-equilibrium case. 

A. Introduction to the Diagrammatic Expansion 

We start by returning to the expansion of 5'(— oo, — oo) in terms of the bare charges in 
Eq. (|l^). We will focus on the expectation value of the Tc ordered product. Let us use a 
slightly different notation, using t's to denote the positions of + charges and s's to denote 
the positions of — charges. Let us take some configuration of charges labeled by tj and Sj, 
with some of them on the top and some on the bottom branch (this way we do not have to 
worry about the superscripts for top and bottom branches, since we can keep track of where 
each charge is by its index). Using this notation, we can write for the Tc bracket: 

where Q = Q\ + Q\ = Q[_ + Qt, and adt, t') = ±1 depending on the ordering of t and t' 
along the Keldysh contour. 

Consider now integer g^s, such that 2g is even and the expression above does not change 
if we take [6 + i(t — t')acit,t')] [6acit,t') + i(t — t')]. The expression for the Tc bracket 
can be simplified with the aid of the following identity which can be proved using partial 



fractions or properties of determinants |20 



iii<j{zi - Zj) n^jK - 

Y{i,j{Zi - Wj) 
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detM{z,w) , (43) 



where M is a matrix defined by 



= ■ (44) 



Zi - Wj 



The presence of the regulators 5 in the expression for the Tc bracket shghtly comphcates 
how we apply the identity to the problem. By naively defining 

dadU, Sj) + i{ti - Sj) ^^^^ 

we would obtain terms in the numerator to order 5 and higher that would not match the 
numerator of the expression for the T,. bracket. This corresponds to a different choice of 
regularization, and we shall return to this point later. The leading term (order 5°), however, 
is exactly the same, and we proceed with the program, writing (detM)^^' for the Tc bracket. 

Now notice that the terms Mij correspond to the interaction between + and — charges, 
and that the expansion of the determinant will be comprised of all ways of combining + to 

— charges in pairs such that each charge only appears once in the expansion. Let us then 
associate a graph for any such dipole combination, as show in Fig. 7. When we raise the 
determinant to the power 2g, the effect is to obtain all different ways to connect the + and 

— charges with hues, such that each charge is connected by exactly 2g lines (see Fig. 7, 
where we illustrate the case of g = 1). 

The graphs so obtained give us a systematic way to account for the contributions to 
Z. The terms in the expansion where all charges are connected to one and only one other 
charge, as in Fig. 8a, are the independent dipole terms. Notice that each line in the graph 
that connects two distant charges roughly corresponds to 1/t. When there are 4 charges, 
the lowest order in 1/t that can be obtained from the expansion comes from taking two 
dipoles and using one line from each to connect it to the other, so that there are 2 lines 
connecting the dipoles (Fig. 8b). In this way we obtain a term which corresponds to 
the leading dipole-dipole interaction. This systematic way to expand Z can be used as the 
formal support for the dipole picture developed in the previous section. 
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B. Equilibrium Case 



To illustrate the power of the formalism described in the previous subsection, we will 
consider the equilibrium case at zero temperature. According to the dipole approximation, 
we expect that the current- current correlation should go as l/t^ for g > 2. In this section 
we will show this is the case for any integer g >2. 

In equilibrium, we no longer need to use the Keldysh contour. Instead, to simplify the 
calculations we will work in Euclidean space, and we will take the field correlator to be 
(O|0(t)0(O)|O) = — ln(t^ + 6^). This differs from the original choice of correlation function 
only in that we have used a different cut-off. The choice given here corresponds to left 
and right movers scattering off each other instead of left movers with left movers, as in the 
original choice. However, for a single impurity, left movers and right movers should really 
be equivalent, so this choice should not make any important difference in the results. More 
importantly, in both cases we choose to regulate the correlator (O|0(t)0(O)|O) consistently, 
no matter where it appears in the expression for P in Eq. (0). It may appear that whenever 
two tunneling events with the same charge interact, we could just ignore the cut-off in the 
numerator of P, since ^e*''"'^'-*-' e*''''^*'^'* ) = {t'^ + is not singular as 5 goes to zero; recall that 
7 = y/g. However, the 5's in this correlator will be multiplied by other correlators that are 
singular as 6 goes to zero, so it turns out that the answer depends on how we regulate the 
numerator. Because we are using the Coulomb gas picture, for now we will choose to keep 
the (5's in the numerator. 

With our choice of regulator, the expression P for the bracket needed to evaluate 
S{—oo, — oo) becomes 



In this equation, the positive charges are at the U and the negative charges are at the Sj. 
Because we are in Euclidean space, we no longer have to use time-ordering when we evaluate 
the integrals over the ti and Sj. 




(46) 
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To simplify the expression for P for any integer g, we will use the same procedure as in 
Ref. |T^. We will write P = AB, where A equals P with the 5's in the numerator set to 
zero, and B is the correction due to the 5's in the numerator of P. Then 

and B is equal to sums over products of 6"^/ (ti — tjY and 6"^/ (sj — SjY, where any one of these 
expressions can occur at most g times in a product. B comes from writing each correlator 
in the numerator as (tj — tj^ [1 + 5^/(tj — tj^] and factoring out the {U — tjY part. 

We can again use the identity in Eq. ( ^31) to simplify the expression for A. If we define 



the matrix M{5) as 



then A is given by 



A = [detMij((5)detMy(-5)]^ . (49) 



As explained in the previous subsection, if we represent each charge by a point and each 
factor of by a directed line, then we obtain all the different ways to connect the 

positive charges to the negative charges so that each charge is connected by exactly 2g lines, 
(half of which are pointing toward the line, and half away from the line). 

We can also give a graphical interpretation of B. Once we have a graph from A, to take 
into account the fact that the numerator is also regulated, we obtain our graphs for P by 
joining any number of pairs of similar charges with the pair of edges l/{ti—tjY or 1/ (sj — s^)^. 
Each of these edges is accompanied by a factor of 5, and any pair of charges can be joined 
by at most g of these pairs of edges. Thus B introduces an interaction between like-charged 
particles. 

In the graphs of A and B, it is important to keep track of the number of vertices, V , 
the number of edges, E, and the number of factors of 5 in the numerator, — /. If we are 
calculating the charge-charge correlation function, and we insert 2N additional charges, 
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then the number of vertices is = 2N + 2. For any connected graph of A, we then have 
E = {2N + 2)g and / = 0. Once we include the effects of B, f is no longer equal to 0, but 
E + f is still given by 



Also, it is important to note that any connected graph is also IPl. This way of describing 
the bracket, P, works similarly in Minkowski space. 

Next, we will evaluate the connected correlation function of (0|e*''"^^*''e~*''"^^*-'|0) for any 
integer g > 1. (The case when g = 1 was considered in Ref. ||19|-) This calculation will also 
work for the correlation functions (0|e^*'''''^*^*'*e^*'''''^*^''^|0), so that these results can be used to 
find the leading dependence on t — s of the current- current correlation functions. 

At the (2A^)*'^ order in perturbation theory, we have 



where A and B depend on t, s, the t^'s and the s^'s. To obtain the connected correlation 
function, we just need to consider the connected graphs in the expression on the right-hand 
side of the above equation. 

To evaluate the integrals, we will perform contour integrals where we complete the con- 
tour in the upper half plane. Thus, for each vertex, tj, we will be evaluating residues for all 
the poles occurring at tj = Sk + iS. (Here, we are using tj and Sk to stand for any type of 
vertex.) We note that in 5, it appears that we will have poles on the real axis. However, 
we know that the original expression for AB does not have any poles on the real axis. This 
means that if we sum over all the graphs for B, these poles cancel, which implies that as 
long as we integrate over the variables in each of these graphs in the same order, we can just 
ignore the poles on the real axis. 

We can describe the process of evaluating residues diagrammatically, as explained in Ref. 

. If the multiplicity of the pole at tj = Sk + iS is equal to one, then there is only one 
edge, Cjk, that joins tj to Sk and represents this kind of pole. In this case, when we evaluate 



E + f = {2N + 2)g. 



(50) 




(51) 
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the residue, we just "collapse" the vertex tj and the edge Cjk- This means we remove the 
vertex tj and edge ejk, and then reconnect all the other edges that were originally connected 
to tj to Sk instead. If the other end-point of any of these edges was also connected to Sk, 
the edge becomes l/{ic5), for some integer c. Otherwise, it remains an edge. Consequently, 
the total number of edges decreases by at least one, and the sum of edges and factors of 6 
in the numerator decreases by exactly one. Also, the graph remains connected and IPI. 

When the multiplicity, m, of a pole is greater than one, then instead of collapsing only 
one edge, we must collapse all the m edges that correspond to the pole. In addition, we must 
take m — 1 derivatives with respect to tj. Each of these derivatives increases the number of 
other legs connected to tj by one, so we obtain m — 1 new legs. Once we have created these 
m — 1 new edges and collapsed both the vertex and the m edges corresponding to the pole, 
we again find that the number of edges decreases by at least one, and the sum, E + f, still 
decreases by exactly one. Again, the graph remains IPI. 

Now we can count the number of edges and factors of 5 that remain after we have done 
all the integrations. The original graph with 2N + 2 vertices has E + f = {2N + 2)g. After 
we integrate over the 2A'" inserted charges, this sum becomes 

^ + /= (2iV + 2)^-2Ar, (52) 

and the only two remaining vertices are t and s. Because the graph must still be IPI, we 
must have at least 2 edges connecting t and s. Since the total number of edges always 
decreases by at least one, we also have 

2<E<{2N + 2)g-2N. (53) 

We will let In = {2N + 2)g — 2N. Lastly, because the final answer must be symmetric in t 
and s, after we sum over all the graphs we can only have even values for E. 

Putting all of this together, we find that the correlation function of (0|re'T^Wr*e-'T^(*)|0) 
must have the form 
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where the a^/s are constants that are determined from integrating the exphcit graphs. 
In order to interpret these results, for the equihbrium case it is helpful to renormalize the 
coupling. We will replace each F and T* with r6^~^ and r*6^^^. This just takes into account 
the self-interaction of the charges and a rescaling of the time variables. 
The correlation function is then 

rr-^ + t (rr-)'-" («.,^ + a.,^ + . . . + . (55) 

This general form is true to all orders in F. Also, note that the derivation of this result did 
not depend on the sign of the charges at t and s, so we will obtain a similar expression for 
two positive charges or two negative charges at t and s. For large times, (or small cut-off S) 
the leading behavior is 

1 00 r2 00 

E (FF*)--^a., + E (rr*)--^a.4. (56) 

This expression appears to go as l/(t — s)^ instead of as the l/(t — s)^ predicted by the 
dipole picture. However, as we shall show shortly, if both the denominator and numerator 
are regulated in the same way, as in Eg. (^Bj) , then = for all A^, so the leading behavior 
does go as S'^/(t — s)^ to all orders in perturbation theory. 

Before showing that ai\i2 = for all A^, we will first use the previous calculation to 
describe a systematic way to determine the leading behavior of each graph. First, we note 
that a final answer of l/(t — s)" corresponds to a graph with n legs joining the vertex t to the 
vertex s. If we remove these n legs, the graph breaks into 2 disconnected components, one 
containing t, and the other containing s. Because the integrations consist only of collapsing 
vertices and edges and also making extra copies of edges, these n legs must have come from 
/ legs in the original graph, where I < n. In addition, because the process of integration 
does not change the connectedness of the graph, when the / legs in the original graph are 
removed, it will break into two disjoint, connected graphs, one containing t and the other 
containing s. An example of this is given in Fig. 9. We also note that since each graph is 
IPI, to break it into two we must remove at least 2 edges. 
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This all implies that the only graphs that can have a leading term of l/(t — s)^ are those 
that are broken into 2 when 2 legs are removed; the only graphs that can have a contribution 
of l/(t — s)^ are those that are broken into 2 when 2, 3 or 4 legs are removed, and, in general, 
only the graphs that can be broken into two connected pieces when 2, 3, . . ., or n legs are 
removed can contribute a term of order — s)". A simple counting argument shows that 
when / legs are removed, the maximum net charge either of the 2 resulting graphs can have 
is l/2g. Because the net charge is always an integer, when / is equal to 2 (and g is greater 
than 1) this means that the net charge must be zero. 

Thus we can classify the graphs according to what their leading behavior is, and we can 
determine which graphs will contribute to any particular term in the expansion in Eg. (p5|) . 
To make contact with the previous subsection, we remark that for the insertion of two 
charges, the only configuration that breaks into two graphs when two lines are removed is 
precisely the one shown in Fig. 8b. 

To arrive at a useful way of estimating graphs (which should also apply in the non- 
equilibrium Minkowski space formalism), we observe that every time we evaluate a residue 
of a pole at = Sj + i6, we are taking tk to be very close to Sj. If we then evaluate an 
Sj = ti + i6 residue, we evaluate Sj close to ti, so in turn that means tj is also close to t/. 
Following through on this observation, we see that for the final result, all the points are 
either evaluated close to t or close to s, and whether it is t or s depends on whether, when 
we remove the n legs, the point is in the graph connected to t or to s. Thus it appears that 
the only contributions to the integral come from all the ways to take some of the vertices 
close to t and the remaining vertices close to s. The exponent of the leading contribution 
will then be determined by the net charge of each of the two resulting subgraphs. This is 
exactly what was done in Section V for the case of two dipoles. 

We now return to calculating the coefficient, aN2, of the l/(t — s)^ part of the charge- 
charge correlator. From the previous discussion, we know that this should come from all 
ways of forming a neutral multipole around t and a neutral multipole around s. As long as 
t — s » 6, we can assume that all the charges in each multipole are much closer to each 
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other than t and s are to each other. We will let to • • • tm-i and sq . . . Sm be the charges close 
to t and tm+i . . . tiv and Sm+2 ■ ■ ■ sjy be the charges close to s. To simplify the notation, in 
most of what follows we will let equal t and Sm+i equal s. Next, we will change variables 
so that ti = Pi + 1, Si = qi + 1 for the charges close to t and tj = pj + s, Sj = qj + s for the 
charges close to s. Then the expression for P becomes 

P = P1P2I', (57) 

where 

p ^ n5=^+M<, {{p^ - vi? + ?>'r - g.)^ + ^'Y .... 

and 



^ _ ni^o Ulm+i {{t-s + p^- p,r + {{t-s + q,- q,Y + 6') 
UT=o nf=^+i {{t-s + p.- q,f + {{t-s + q,- p,f + • 



(60) 



Pi and P2 just look like the original integral, but for a smaller number of charges, so they 
contain the intra-multipole interactions. The expression for / contains all the interactions 
between the two different multipoles. In the numerator, the positive and negative charges 
of the first multipole interact with charges of the same sign in the second multipole, and 
in the denominator the charges of the first multipole interact with the charges of opposite 
sign in the second multipole. Because the multipoles are both neutral, and because every 
factor in Eqn. ( |60|) depends on t — s, both the numerator and denominator have the same 
number of factors oit — s. Once we divide through by t — s, similar counting tells us that the 
number of times Pi/{t — s), qi/{t — s) and — s)^ each appear in the numerator equals 
the number of times each of these appear in the denominator. If we expand / out for large 
t — s and count all the terms that contribute to order — s)^, we find 

^ = 1 + 7ri-T^2 Y^iPiPi + m - Pilj - Pjli)^ (61) 
^) i,j 
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where pi and qt run over all the charges in the first multipole and pj and qj run over all the 
charges in the second multipole. Then 

= '^ + rrhvi ^^p^p^ + ^^^^^ ~ p^'i^ ~ Pi'i^^- (^2) 

-^J i,j 

The important feature of the — s)^ part is that it is odd under changing the signs of the 
coordinates of all the charges in only one multipole. Meanwhile, P1P2 is even under such a 
sign change, so once we integrate over all the coordinates, the l/{t — sY part vanishes and we 
are left only with the — s)^ (which should vanish once we sum over all configurations of 
the charges) and the — sY parts. Thus, the coefficients, a2N should vanish to all orders 
in perturbation theory and the charge-charge correlation functions, (0|re*''''^''*^r*e~'''''^''*^|0), 
should go as 04^^/ {t — s)^, for some constant 04. It is considerably more difficult to evaluate 
this constant. 

One final remark is that if we had regulated only the denominator, then the previous 
argument would not have gone through: the — s)^'s from the denominator would no 

longer be canceled by the — s)^'s from the numerator, so that a2N would be non-zero. 
In this case, the correlation functions instead would go as l/(t — s)^. 

C. Implications for the Non-equilibrium Case 

Even for the non-equilibrium case, we can use our analysis of the graphs in the preceding 
subsection to guide us in determining which graphs should give the leading contributions 
to the current-current correlation functions. To calculate the singularity at = 0, we can 
use the same neutral multipole expansion as in the end of the previous section. The only 
changes to Eqs. ( p8| , |59| , and ^0|) for the intra-multipole and inter-multipole interactions 
are that we must now use the non-equilibrium regulators which depend on the ac(tj,Sj)'s. 
Also, Eq. (IID for the muhipole Pi will now be muhiplied by n^o e*"^"^' Tl^o e"'"^"^" and 
Eq. dl) for P2 will be multiplied by Ililm+i e*"""^' nJLm+i e"''^''^'- Consequently, P1P2 no 
longer remains unchanged when all the signs of the vertices in one multipole are reversed. 
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Therefore, according to Eq. (|62D the contribution to the current-current correlation function 
when one multipole is close to vertex t and the other is close to vertex s goes as 

fc=o 1=0 ij 

where pi and qi are in the first multipole and pj and qj are in the second multipole. Also, 
we only take the connected graphs in the multipoles Pi and P2- Thus, to all orders in F, 
the correlator goes as l/(t — s)^ + O {l/\t — s\^). This means that, at low frequency, the 
noise spectrum should have a singularity that goes as \u!\ at every order in F. Here we 
are assuming that for (7 > 1 the neutral multipoles are all bound, just as they are in the 
equilibrium case. 

In the non-equilibrium case, we also expect singularities at = ±uJo and possibly also at 
uj = nujQ for other integer values of n. To find the leading behavior at these singularities we 
use the fact that the expression for P in Eq. (^) can be expressed as a product AB. As in 
the preceding subsection, A is a determinant, and B contains the corrections that naively go 
as (1 + 0(5)). For the non-equilibrium case, A was defined at the beginning of this section 
as det Mij where 

Mii = 1-^ r^— 7 r. (64) 

dac{ti,Sj) + i{ti - Sj) 

The graphs for the A defined here are identical to those in the previous subsection, except 
for the choice of regulator. This means that all of our previous counting arguments should 
apply. However, the form of B is now much more complicated than before, so it is not clear 
whether it modifies the counting in the same simple way as before. Because the expression 
for P in Eq. (^) and the expression for A both contain the information about which branch 
each charge is on, and since the only difference between the two expressions is the choice of 
regulator, for convenience we will choose to work with A = det Mij instead of with P. (In 
the equilibrium case, we have seen that picking a different regulator does not change the 
types of terms that can appear in the final answer; it just changes the value of the coefficient 
in front of each term, possibly setting some to zero. In case of a discrepancy, the choice of 
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regulator should reflect the physics at hand, so it is useful to keep in mind that in P the 
interactions in the Coulomb gas are regulated and in A the fermion-like propagators in the 
matrix M are regulated.) 

For detMjj, our counting and classification of graphs proceeds as before. This implies 
that if we can break the graph into two connected multipoles with charge Q and —Q, 
respectively, then the graph will give a leading contribution of 

^iQuJot^—iQuJos 

as long as all charges within a multipole are close to one another. This will give the singu- 
larity \uj ± QcuoP*^^"^- For example, the graph in Fig. 9b will give a contribution as in Eq. 
(|65|) with Q = 1 and g = 2. Without performing the integral, we cannot determine whether 
ag (which can depend on 6 and Uq) is non-zero. However, from this line of reasoning, we 
can conclude that the rr*|co' ± uqI"^^'^ singularity should only receive corrections that go 
at least as Itu ± ujq\'^^~^ at all higher orders in FF*. Similarly, at higher multiples of ujq we 
expect the singularities to be even smoother because they go at least as Icj ± QuJol"^^^'^. 

As a check on these calculations, we note that we can apply the same analysis of the 
graphs and similar counting arguments even at g = 1. In this case, every connected graph 
is just a simple polygon with alternating charges at the vertices. It is straightforward to 
see that when any such graph is divided into two disjoint, connected parts, each part can 
only have a total charge of or ±1, and exactly two lines must be cut. Therefore, the only 
singularities we can obtain are \uj\ and Itu ± cuqI, with no higher order corrections. These 
results agree with those in Ref. and give strong evidence that our method of analyzing 



the graphs works even for the non-equilibrium case. 

VII. CONCLUSION 

In this work we defined a framework for the study of equilibrium and non-equilibrium 
noise in 1-D Luttinger liquids. The interactions give rise to correlations that are manifest in 
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the noise spectrum. The correlations are responsible both for algebraic singularities in the 
noise power spectrum and for the nonlinear dependence of the strength of such singularities 
on either the applied voltage between the terminals of the 1-D system or the temperature. 
The information carried by both the form of the singularities and their strength can help us 
identify Luttinger liquid states in experiments. 

The picture of the tunneling in terms of the Coulomb gas (and its dipole gas interpre- 
tation) is attractive because it gives us an intuitive way to think about the tunneling in 
the Keldysh formalism. This picture provides a unified description of the low and high 
frequency noise: correlations between different dipoles define the structure of the noise near 
zero frequency, whereas correlations between the two charges within the dipole should con- 
tribute to the noise near the "Josephson" frequency cuj — e*V/h. Using formal diagrammatic 
techniques we have justified this interpretation, and, for integer g, we have obtained exact 
answers for the form of the singularity in the equilibrium case. 

One particulary striking result we obtained is that the form of the leading singularity 
at zero frequency (oc |a;|) is the same for strongly correlated Luttinger liquids as well as 
for non-interacting systems. The effects of correlations in the case of low-frequency noise is 
present only in the strength of the singularity, with a strong non-linear dependence on the 
applied voltage that is proportional to V^(s-i) 

Although our Coulomb gas picture and the accompanying formalism has enabled us 
to calculate the form of the singularities to all orders in perturbation theory, beyond the 
order |r|*^ it is too cumbersome to find the strength (i.e. the coefficient in front) of these 
singularities. We would also like to point out that the structure of the noise far away from 
the frequencies ncuj, at higher orders in perturbation theory, is unknown; the information 
we are able to obtain is limited solely to frequencies near the singular points. One exception 
is the exactly solvable case of g — 1, where we find that the noise spectrum must have the 
form a + b\u!\ + c\u! ^ u)j\, where a, b and c can be calculated from the non-equilibrium 
voltage and the transmission coefficient. Thus in this case we recover the results for non- 
interacting electrons. Indeed, the framework we presented can be used with = 1 for 
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studying coherence effects which appear in the noise for non-interacting electrons and are 
due to the Pauh principle, because the statistics enter in the formulation we use through 
the bosonization. 

There are two points in this work that need further exploration. The first is the apparent 
fine point of better understanding the role of the short distance cutoff in our calculations. 
We need either to determine whether the non-equilibrium noise is sensitive to our choice of 
regulator or else to show that our choice of regulating the fermion-like propagators instead 
of the Coulomb gas is the physical one. The second, and more important, question is to 
understand non-pertubative effects. For example, one expects that the position of finite 
voltage singularities should depend on F. In the case of tunneling between edge states, 
when we increase the current, the frequency should shift from to as we go from the 
configuration in Fig. Ic, where the electrons are tunneling, to the one in Fig. lb, where 
the quasi-particles are tunneling. This is not reflected in our perturbative calculations. 
However, we have evidence that within our Coulomb gas picture this shift can be explained 
by non-perturbative effects, and we hope to address this issue in a future paper. 
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FIGURE CAPTIONS 



Figure 1. Schematic drawing of the geometries for tunnehng in 1-D Luttinger Uquids. A 
channel connected to two reservoirs is shown in (a), with a potential barrier or weak link in 
the middle. The geometries for tunnehng between edge states are shown in (b) and (c) . By 
adjusting the gate voltage Vq one can obtain either a simply connected QH droplet (b), or 
two disconnected QH droplets (c). For the geometry in (b) both electrons and quasiparticles 
(carrying fractional charge) can tunnel from one edge to the other, whereas for the geometry 
in (c) only electrons can tunnel. The tunneling current It depends on the applied voltage 
between the right and left edges. 

Figure 2. An insertion of an operator e'*"*''"^^*^ correspond to the insertion of a charge 
+ on the contour at time t. Similarly, an insertion of an operator e"*''"^^*-* correspond to 
an insertion of a charge — at time t. The time t is ordered along the contour shown, 
and there is a distinction between charges placed on the top and bottom branches. For 
illustration, in the example shown we have for the number of + and — charges in the t and 
b branches Q\. — 3, = 2, = 2 and Qt. — 3. Only terms that have zero total charge 
Q — Q\ + <5+ — Q^- — Q^- can give a non-zero contribution to Z. 

Figure 3. The apphed voltage V between the terminals or edges creates an unbalance of 
charge between the top and bottom branches. Since + and — charges correspond respectively 
to tunneling from i? — > L and L — > i?, an excess of charge in the top branch correspond to 
net tunneling in one direction. 

Figure 4. The charges that form the Coulomb gas can form a dipole phase. In this phase, 
the expression for Z can be recast as a sum over dipole strengths and positions, instead of 
summing over the locations of the + and — charges. 
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Figure 5. The four types of dipole, classified according to the position of the + and — 
charges comprising it. In the c_|_ dipole the + charge is on the top branch and the — charge 
is on the bottom. In the c_ the — charge is on the top and the + is on the bottom. In 
the t dipole both charges are on the top branch, and in the b dipole both charges are on 
the bottom branch. Notice that only the c± dipoles contribute to a net current, as they 
create an unbalance of charge between the top and bottom branches. The t and b dipoles 
contribute to the noise, but not to the current. 

Figure 6. Two dipoles will interact because of the relative position between the charges 
that comprise them. The figure shows two dipoles with center of mass positions ti and t2 
and strengths pi and p2- 

Figure 7. The expression for the correlation between many charges can be expressed as 
a power of the determinant of a matrix M. The matrix element Mij can be represented 
diagrammatically as a line connecting a + charge at position ti to a — charge at position Sj, 
as shown in (a). The determinant contains different ways of pairing the charges (b). Finally, 
when raising the determinant to the power 2g (done in this figure for = 1), we generate 
different ways of connecting the charges such that exactly 2g lines leave each + charge and 
exactly 2g lines arrive at each — charge, as shown in (c). 

Figure 8. The graph corresponding to independent dipoles is shown in (a), with all lines 
leaving the + charge arriving at the — charge (here we use = 2 for illustration) . One of the 
graphs contributing to the dipole-dipole correlation is shown in (b). Each leg connecting the 
two dipoles contributes to order 1/t, so that the dipole-dipole correlation is of order 

Figure 9. A sample graph with g = 2 that gives a contribution of l/{t — s)^ after it is 
integrated. The final graph with 4 legs is shown in (a). It is obtained by integrating over 
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the vertices ti, t2, si and S2 in the graph shown in (b). The 4 final legs in the final graph 
come from the 4 bold-faced legs. In (c), the two disjoint graphs (or multipoles) resulting 
from removing the 4 bold-faced legs are shown. 
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